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To identify the possible mechanism of coercivity (H c ) degradation of Nd-Fe-B sintered magnets, we study 
the roles of the exchange field acting on the 4f electrons in Nd ions and theoretically investigate how the 
variation of the exchange field affects the values of the magnetic anisotropy constants K\ and K 2 . We find 
that, with decreasing exchange field strength, both values decrease as a result of the lower asphericity of 
the 4f electron cloud, indicating that the local anisotropy constants might become small around the grain 
boundaries where the exchange fields are decreased owing to the smaller coordination number. 


Nd-Fe-B sintered magnets 1_3 ' have the largest maximum energy product among the current mag¬ 
nets and have been widely used for magnetic devices such as voice coil motors in magnetic recording 
systems. Recently, because of the rapidly growing interest in electric vehicles, much effort has been 
made to suppress the degradation of the coercivity ( H c ) of Nd-Fe-B magnets. However, from an in¬ 
dustrial viewpoint, reduction in the usage of Dy is strongly desired, because Dy is a rare metal and 
the magnetization of the Nd-Fe-B magnets decreases by substituting Dy with Nd owing to the antipar¬ 
allel coupling between Dy and Fe moments. Realizing Dy-free high-performance Nd-Fe-B magnets 
requires a further increase of H c in the Nd-Fe-B system by microstructure optimization, 4 ^ 11 ' and there¬ 
fore, establishing the microscopic foundation for the coercivity mechanism is desired. From a theoret¬ 
ical viewpoint, many works 12 ' l6 ' have focused on the change of magnetic anisotropy constants around 
the grain boundary surfaces as a result of the stresses, defects, and change of spatial symmetry. In 
addition, micromagnetic model calculations have shown that the surface c-plane anisotropy can dras¬ 
tically decrease the coercivity. 17 ' For these reasons, evaluation of the local anisotropy constants around 
the grain boundaries and determining their temperature dependence are important to investigate the 
degradation of the coercivity. 

With regard to the magnetic anisotropy of rare earth (RE) transition metal compounds, it is be¬ 
lieved that the 4f electrons in RE ions are responsible for the main part of the magnetic anisotropy and 
that the crystalline electric field (CEF) acting on the 4f electrons dominates this property. 18 ' Under the 
assumption that the exchange field on the 4f electrons is strong enough, by using the CEF parameter 
A?,', the leading anisotropy constant K\ can be approximately described by K\ = -3J(J- l)a(r 2 )A®NR, 
where a is the Stevens factor, J is the total angular momentum, (r 2 ) is the average of r 2 over the radial 
wave function of the 4f electrons, and Nr is the density of RE ions. Note that the CEF parameter 
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is easily affected by circumstances especially around the grain boundaries. Actually, Moriya et al. 15) 
predicted by using a first-principles calculation that the CEF parameter A® exhibits a negative value at 
the (001) surface of Nd 2 Fei 4 B . According to Ref. 17, this may lead to a drastic degradation of H c . 

Worth mentioning here is that the exchange field acting on the 4f electrons from surrounding Fe 
spins can also be changed as a result of the situation surrounding the RE ion. It is natural to consider 
that the exchange fields arc weak around the grain boundaries compared to those in the bulk. In this 
case, the above expression for K\ no longer applies because the precondition that the exchange field is 
strong enough no longer holds. Especially, the decrease of the exchange field possibly has a significant 
influence on the temperature dependence of the anisotropy constants because of thermal fluctuations 
of the 4f moments. Based on this viewpoint, we focus on the roles of the exchange field on the 4f- 
related anisotropy constants, and we theoretically investigate how the magnetic anisotropy is affected 
by a variation of exchange fields and the CEF acting on the 4f electrons. To this aim, the temperature- 
dependent anisotropy constants K\{T ) and Ki{T) arc calculated based on crystal field theory. 

The calculation method we use here is based on conventional crystal field theory for the 4f elec¬ 
tronic system in the RE ions. The total Hamiltonian for 4f electrons is given as 77 = 7T K)n + 'Hew + 
77 ma g, where 7T K)n describes the intra-atomic interactions in an RE ion, TTcef denotes the CEF pari 
representing the electrostatic field acting on the 4f electrons from the surrounding charge, and 77 ma „ 
denotes the effective exchange interactions between the 4f and Fe moments. When Nd ions arc in¬ 
volved, only the ground J multiplet can be considered as the 4f states because 4f electrons of the Nd 
ion have large spin-orbit coupling compared to the CEF splitting. In this case, 7 /cef can be written 
as ‘Hr.F.p = X/,m B'J l d’", where B'" = f){{r l )A'j 1 is the CEF coefficient and O” is the Stevens operator 
expressed by the multinomial of the orbital angular - momentum operators. 9j is the Stevens factor a, 
(i, or y for / = 2, 4, or 6, respectively, and A™ are the CEF parameters. In this work, we take into 
account only B®, B% and B®, for simplicity. Within the ground J multiplet, 7i ma „ can be expressed as 
"Knag = 2(gj - 1)./ H ex , where gj is the Lande factor and H ex is the effective exchange field reflecting 
the exchange interaction with surrounding Fe spins. In principle, H ex depends on temperature through 
a self-consistent equation for a molecular - field approximation in the full system of Nd 2 Fei 4 B. How¬ 
ever, in the present work, we assume that H ex is proportional to a molecular - field of Fe spins (Sp e ) 
and that the molecular - field is described by the Heisenberg model for the Fe-spin system. Solving 
the self-consistent equation for (Sp e ) under the condition that the spin of the Fe ion is 1 and that the 
Curie temperature T c is 583 K, we have confirmed that (Sp e ) reproduces well the temperature depen¬ 
dence of the magnetization of Y 2 Fei 4 B. 19) Thus we have introduced temperature dependency in H ex 
by H ex (T ) = Hex (Spe), where H ex is the strength parameter of the effective exchange field. Here, the 
direction of total magnetization is assumed to coincide with that of H ex , because the total magnetiza¬ 
tion is mostly governed by the Fe moments even when the direction of 4f moments differs from that 
of the Fe moments. 

By defining 9 as the angle between H ex (T ) (total magnetization) and the c axis of the crystal 
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lattice, the free energy of the 4f electronic system can be written as 


F(0,T) = -k B T InTrexp 


< H{0, T) 
k B T 


where the eigenvalues arc calculated by diagonalizing 'H(0, T) within the J = 9/2 subspace (10 x 10 
matrix), and F(0, T ) is numerically obtained. 

The 4f-related anisotropy constants K\ and K 2 arc coefficients of sin 2 6 and sin 4 0 in F(6, T). For 
small 0, F(0, T ) can be expanded as 


0 ~ F 4 

F(0, T ) = F(0, T) + —F a \0, T ) + —F ( 4 ) (0, T) + 


= F( 0, T) + K l (T)0 2 + 


-^Ki(T) + K 2 (T) 


06 

—F (6) (0, T) + ■ ■ ■ 
6 ! 

U 4 + • • •. 


Therefore, K\ and K 2 can be expressed by 


Ki(T) = ^F ( 2 ) (0,F), 

K 2 (T)= X -K { + j [ F (A \0,T), 

where we calculate F {2 \ 0, T) and F ( 4 , (0, T) from finite differences. 

To perform a quantitative analysis on the magnetic anisotropy constants of Nd 2 Fei 4 B, we start by 
reproducing the experimental results of K\ and K 2 of bulk Nd 2 Fc i jB by using the present approach. 
Figure 1 shows the calculated results obtained for K\(T) and K 2 (T ) using A { \ = 450 K/ao 2 , A 6 = 
-45 K/ao 4 , A® = -0.1 K/ao 6 , and H ex = 364 K, 20 ’ together with the experimental data . 21 - 1 Here, ao 
is the Bohr radius. In calculating K\{T), we add the Fe sublattice contribution K ] { C (T) deduced from 
experimental data for Y 2 FC 14 B. I9) This has an almost constant value of MJ/m 3 below 300 K and 
decreases to zero when the temperature approaches T c . The difference between two Nd sublattices is 
not taken into account here. The inset in Fig. 1 shows the tilting angle of the magnetization vector 
calculated from (>F( ^ T) = 0. Including this spin reorientation behavior, one can see agreement between 
calculated data and experimental ones to a certain degree. 

It should be noted here that the experimental values of K\(T ) and K 2 (T ) arc not obtained by the 
direct measurement but are usually deduced from magnetization curves or torque curves under finite 
applied fields. Furthermore, in our numerical analysis, we made some approximations to express K\(T) 
and K 2 {T). In this sense, we believe that the CEF and H ex parameters reported by Yamada et al. arc 
more reliable, because in their analysis the Zeeman term is considered in addition to the Hamiltonian to 
ascertain the magnetization curves and to make a direct comparison with the measured ones. Reflecting 
this difference, the CEF and // ex parameters obtained by Yamada et al. arc considerably different from 
ours. To examine the availability of the present method and reliability of the parameters obtained 
above, we also calculated the K\{T) and K 2 (T ) curves using the parameters given by Yamada et al. l8) 
Except for the behavior of K 2 (T ) below 70 K, we confirmed that the calculated data agree with both 
the experimental and our calculated results. These results lead us to believe that the present analysis 
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Fig. 1 . Calculated anisotropy constants K\ and Ko (solid lines) as a function of temperature obtained using 
A° = 450 K/ao 2 , A® = -45 K/ao 4 , A® = -0.1 K/ao 6 , and H ex = 364 K. The closed and open circles are the 
measured values 21 ' of K\ and K 2 , respectively. The inset shows the tilting angles of the magnetization from the 
c axis as a function of temperature together with the experimental results (closed circles). 21 ' 
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Fig. 2. Calculated anisotropy constants K\ and /C (solid lines) as a function of temperature obtained using 
the same parameters as in Fig. 1 (dashed lines) but with the half value of // ex . 


and the parameters obtained here can provide an appropriate characterization for K\(T) and K^iT) at 
the quantitative level at least above 100 K. 

Next we proceed to examine how K\(T) and Ki{T) arc affected by the variation of amplitude of 
H ex (T). Figure 2 shows the calculated K\(T) and Ki{T) using the half value of H ex with the other 
parameters unchanged. The results in Fig. 1 arc also shown for comparison. One can see a steep 
reduction in both K\(T) and Ki{T) near - a temperature of 70 K. This in turn leads to the result that, for 
temperatures above 200 K, both K\ and Ki values become smaller than those in Fig. 1. To see the H ex 
dependence of K\ and Ki values, we plot in Fig. 3 these values as a function of H ex at T = 300 K. 
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Fig. 3. H ex dependence of the anisotropy constants K\ and K 2 at T - 300 K. The CEF parameters are the 
same as in Fig. 1. 

In these calculations, we neglect the Fe contribution K^ e (T) to focus on only the effects of H ex on the 
4f-related anisotropy constants. One finds that K 2 exhibits a monotonic increase with H ex whereas K\ 
has a peak at around H ex = 300 K and goes into the negative region above H ex = 800 K. The negative 
value of K\ for the high-// ex region is naturally understood by considering the fact that K\ starts from 
a negative value at T = 0; that is, the high limit of H ex effectively corresponds to the low-temperature 
limit. Further attention should be paid to the low -H ex region below H ex = 300 K, where both K\ and 
Ki values degrade with decreasing H ex . This is because the 4f electron cloud approaches a spherical 
one, in accordance with the decrease of H ex , resulting in an insensitivity to the CEF and a decrease 
of the anisotropy energy. From this behavior, one should recognize that the temperature dependence 
of the anisotropy constants K\ and K 2 comes not only from the thermal fluctuation of the 4f moments 
but also from the strength of H ex (T ) through its temperature dependence as (Spe)- This can explain 
the vanishing of K\, K 2 , and hence H c at the Curie temperature. Further, even at lower temperature, 
the anisotropy constants may become small around the grain boundaries where the exchange fields 
become small owing to the smaller coordination number or to some other factors such as stresses and 
lattice defects. For these reasons, at around the operating temperature of electric vehicle motors (500 
K), the decrease in the amplitude of H ex in addition to the reduction of (5p e ) owing to its temperature 
dependence may significantly influence magnet performance. 

Finally, we investigate the effects of variation of the CEF parameters on K\(T) and K 2 (T). As 
pointed out in our previous work, 151 there may be a possibility that A 11 exhibits a negative value at 
the grain surfaces, depending on whether their characteristics differ from those of the bulk. Figure 
4 shows the behaviors of K\(T) and K 2 (T) using the parameters A® = -820 K/ao 2 , A^j = A° = 0, 
and H ex = 364 K. Clearly seen is that K\(T) exhibits a negative value for any temperature and that 
K 2 (T) vanishes at temperatures above 200 K, which implies that the system favors planar anisotropy 
above 200 K. Thus the prediction of K < 0 when A® < 0 reported in the previous study is justified in 
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Fig. 4. Calculated anisotropy constants K\ and K 2 as a function of temperature obtained using A® = 
-820 K/a 0 2 , A° = A° = 0. and H cx = 364 K. 


the room-temperature region, from which Mitsumata et al. I7) suggested that H c drastically degrades 
owing to the negative K at the surface. 

To conclude this paper, we summarize our study as follows. To identify the possible mechanism 
of coercivity (H c ) degradation of Nd-Fe-B sintered magnets, we have investigated how the magnetic 
anisotropy is affected by a variation of exchange fields and the CEF acting on the 4f electrons. To 
this aim, the temperature-dependent anisotropy constants K\ (T) and Ki{T) were calculated based 
on crystal field theory. Using a certain set of CEF parameters and exchange field strength, we can 
obtain K\ (T) and KjiT) curves that reproduce well the experimental data for bulk Nd 2 Fei 4 B. It was 
found that, with decreasing exchange field strength, both values decrease as a result of the lower 
asphericity of the 4f electron cloud. This feature makes K\{T), Ki{T), and hence H c vanish at the 
Curie temperature, where the exchange field strength tends to zero. Even at a lower temperature, the 
anisotropy constants may become small around the grain boundaries where the exchange fields become 
small owing to the smaller coordination number. For these reasons, at around the operating temperature 
of electric vehicle motors (500 K), the decrease in the amplitude of H ex in addition to the reduction 
of (Sp e ) owing to its temperature dependence may significantly influence magnet performance. It is 
also confirmed that the negative value of A®, the leading CEF parameter, results in a planar anisotropy 
at room temperature. Thus the prediction of K < 0 when A, < 0 reported in the previous study is 
justified at the room-temperature region. 
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